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O ■ Abstract 

<D ' 

We review two inter-related topics. First, we consider the behaviour of "soft" scattering 
observables, such as crtot> da e \/dt, dasD / dtdM 2 , particle multiplicities etc., at high-energy 
proton- (anti)proton colliders. We emphasize the sizeable effects of absorption on high- 
energy 'soft' processes, and, hence, the necessity to include multi-Pomeron-Pomeron in- 

00 i 

teractions in the usual multi-channel eikonal description. We describe a multi-component 
model which has been tuned to the available data for soft processes in the CERN-ISR 
to Tevatron energy range, and which therefore allows predictions to be made for 'soft' 
observables at the LHC. The second topic concerns the calculation of the rate of exclusive 
processes of the form pp — > p + A + p at high energy colliders, where A is a heavy system. 
In particular, we discuss the survival probability of the rapidity gaps (denoted by the + 
signs) to both eikonal and enhanced soft rescattering effects. At the LHC energy, the most 
topical case is when A is a Higgs boson. At the Tevatron, measurements have been made 
for the exclusive diffractive production of various systems: A being either 77, dijet, or 
a Xc meson. We compare the observed rates with the expectations. Finally, we describe 
how predictions for exclusive processes may be checked in the early runs of the LHC. 



1 Motivation for revisiting soft interactions 

The description of the high energy behaviour of "soft" scattering observables such as <7 to t, 
da e \/dt, dasn / dtdM 2 , particle multiplicities etc. at hadron colliders predated QCD. It is some- 
times regarded as the Dark Age of strong interactions. However, this is not totally fair. There 

1 Based on a talk by Alan Martin at the Cracow Epiphany Conference on "Hadron interactions at the dawn 
of the LHC, 5-7 January 2009; dedicated to the memory of Jan Kwiecinski". 



was a successful description of these processes in terms of the exchange of Regge trajectories 
linked to particle states in the crossed channels [1]. The dominant exchange at high energy 
is the Pomeron, and we have Gribov's Reggeon calculus [2] to account for the multi-Pomeron 
contributions. However the available data did not reach high enough energy to distinguish 
between the different scenarios [3, 4] for the high-energy behaviour of the interaction amplitude 
[5, 6]. 

With the advent of the LHC, it is valuable to revisit soft interactions 2 because of 

• the intrinsic interest in obtaining a reliable, self-consistent model for soft interactions, 
which may be illuminated by data from the LHC; 

• the need to use the model for predictions of the gross features of soft interactions; in 
particular, to understand the structure of the underlying events at the LHC; 

• the advantages of using central exclusive production, pp — > p + H + p, to study the 
properties of the Higgs sector in an especially clean environment at the LHC [9] - [14]. 
For example, the Higgs mass(es) can be measured with very good accuracy (AM# ~ 
1 — 2 GeV) by the missing-mass method by detecting the outgoing very forward protons. 
However, the exclusive cross sections are strongly suppressed by the small survival factor, 
5 2 <C 1, of the rapidity gaps. Thus we need a reliable model of soft interactions to 
evaluate the corresponding value of S 2 [15, 16]. 

2 Description of soft interactions 

We start with Fig. 1, which shows sketches of the optical theorems for the total cross section 
and for high-mass single diffractive dissociation. The high-energy expressions given for the cross 
sections arise from the Regge behaviour of the bare amplitudes. However there are screening 
(that is, absorptive) corrections, which will suppress the cross sections. For example, the value 
of the triple-Pomeron coupling, g$p, originally obtained [17] using the bare formula, will, on 
taking account of the absorptive corrections [18] , be increased. 

2.1 Single channel eikonal — elastic scattering 

The total and elastic cross sections are usually described in terms of an eikonal model. At 
high energy the position of the fast particle in the impact parameter, b, plane is, to a good 
approximation, frozen during the interaction. Thus the value of b is related to the orbital 

2 Two papers resulted from Jan Kwiecinski's first extended visit to Durham in 1990. These were entitled 
"Partons at small x" [7] and "Semihard QCD expectations for pp (or pp) scattering at CERN, the Tevatron 
and the SSC"[8]. It is interesting to observe that the first rapidly developed into a rich and very enjoyable 
collaborative programme, with impact on HERA, whereas the second was almost 20 years premature. 
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Figure 1: Optical theorems for the total cross section and for high- mass diffractive dissociation. 
The high-energy Regge expressions shown for the bare amplitudes have sizeable absorptive 
corrections. 



angular momentum I = by/s/2 of the incoming hadron. Elastic unitarity, therefore, takes the 
form 3 

2ImT el (s,6) = \T cl ( Sl b)\ 2 + G incl (s,b), (1) 

from which it follows that 

^ = 2ImT el M) = 2(l-e^/ 2 ) (2) 

d(T e l 
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where Q(s, b) > is called the opacity (optical density) or eikonal 4 . It is the Fourier transform 
of the two-particle (s-channel) irreducible amplitude, A(s,qt). That is 



|T el (,,6)| 2 = (l-e-^ 2 ) 2 , (3) 
2ImT el (s,6) - \T cl (s,b)\ 2 = 1 - e~ n , (4) 



5 

fi(«,6) = ^ Jd 2 q t A(s,q t )e^ h , (5) 

where q\ = —t, and where the amplitude is normalized by the relation cr tot (s) = 2ImT e i(s, t — 0). 
^From (4), we see that exp(— ft(s, b)) is the probability that no inelastic scattering occurs at 
impact parameter b. In the framework of the eikonal model, the elastic amplitude, is obtained 
by the sum of Regge-exchange diagrams, as shown in Fig. 2(a), which is equivalent to the 
iteration of the elastic unitarity equation, (1). 

At high energies the ratio ReT el /ImT el is small and can be evaluated via the dispersion 
relation. So the imaginary part of Q is usually neglected. 



3 Ginei accounts for the presence of inelastic channels. 
4 Sometimes 0/2 is called the eikonal. 



3 We use the bold face symbols qi and b to denote vectors in the transverse plane. 
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Figure 2: (a) The single-channel eikonal description of elastic scattering; (b) the multichannel 
eikonal formula which allows for low-mass proton dissociations in terms of diffractive eigenstates 
\<pi), \4>kj] and (c) the inclusion of the multi-Pomeron-Pomeron diagrams which allow for high- 
mass dissociation. 



2.2 Inclusion of diffractive dissociation 

So much for elastic diffraction. Now we turn to inelastic diffraction, which is a consequence 
of the internal structure of hadrons. Besides the pure elastic two-particle intermediate states 
shown in Fig. 2(a), there is the possibility of proton excitation, p — > N*. At high energies, where 
the lifetime of the fluctuations of the fast proton is large, r ~ E/m 2 , the corresponding Fock 
states can be considered as 'frozen'. Each constituent of the proton can undergo scattering and 
thus destroy the coherence of the fluctuations. As a consequence, the outgoing superposition 
of states will be different from the incident particle, so we will have inelastic, as well as elastic, 
diffraction. 

To discuss inelastic diffraction, it is convenient to follow Good and Walker [19], and to 
introduce states \<f>k) which diagonalize the T matrix. Such, so-called diffractive, eigenstates 
only undergo elastic scattering. To account for the internal structure of the hadronic states, we, 
therefore, have to enlarge the set of intermediate states, from just the single elastic channel, and 
to introduce a multichannel eikonal. The situation is pictured in Fig. 2(b). It is straightforward 
to show that 

^ 2 (D md ^ = n > (6) 

where the brackets of (T) mean that we take the average of ImT(s, b) over the initial probability 
distributions of diffractive eigenstates of the 'beam' and 'target' protons. The cross section for 
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the single diffractive dissociation of the 'beam' proton, 



= <(T) 2 > - (T) 2 , (7) 

is given by the statistical dispersion in the absorption probabilities of the diffractive eigenstates. 
Here, the average is only taken over the diffractive components into which the 'beam' proton 
dissociates. 

At first sight, it appears that if we were to enlarge the number of eigenstates \<f>i), then we 
may include even high-mass proton dissociation. However here we face the problem of double 
counting when partons originating from dissociation of the beam and 'target' initial protons 
overlap in rapidities. For this reason high-mass dissociation is usually described by "enhanced" 
multi-Pomeron diagrams. The first and simplest is the triple- Pomeron graph, shown in Fig. 1, 
and again in Fig. 2(c). 



2.3 The importance of absorptive effects 

Already at Tevatron energies the absorptive correction to the elastic amplitude, due to elastic 
eikonal rescattering, is not negligible; it is about —20% in comparison with the simple one 
Pomeron exchange. After accounting for low- mass proton excitations (that is iV*'s in the 
intermediate states) the correction becomes twice larger (that is, about —40%). Indeed, the 
possibility of proton excitation means that we have to include additional inelastic channels 
which were not accounted for in the irreducible amplitude A of (5). This enlarges the probability 
of absorption for the elastic channel, that is the effective opacity f2. 

Next, when we account for the screening of high-mass diffractive dissociation, dasu/dM 2 , 
there is an extra factor of 2 coming from the AGK cutting rules [20]. Thus, the absorptive 
effects in the triple-Regge domain are expected to be quite large. The previous triple-Regge 
analyses (see, for example, [17]) did not allow for absorptive corrections and the resulting triple- 
Regge couplings must be regarded, not as bare vertices, but as effective couplings embodying 
the absorptive effects [21]. Since the inelastic cross section (and, therefore, the absorptive 
corrections) expected at the LHC are more than twice as large as that observed at fixed-target 
and CERN-ISR energies, the old triple-Regge vertices cannot be used to predict the diffractive 
cross sections at the LHC. 

Thus, it was necessary to perform a new triple-Regge analysis that includes the absorptive 
effects explicitly. Such an analysis has recently been performed [18] in which the fixed-target 
FNAL, CERN-ISR and Tevatron data, that are available in the triple-Regge region, are fitted in 
terms of l PPP\ l PPR\ l RRP\ l RRR and nnP contributions 6 . To account for the absorptive 

6 The analysis of [18] had the limited objective of estimating the triple-Reggeon couplings; in particular the 
triple-Pomeron coupling g^p = gppp- It should not be used as a model to describe soft high energy interactions 
in the whole rapidity interval. Section 2.4 describes such a model. 
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corrections a two-channel (Good- Walker) eikonal model was used, which describes well the 
total, a tot , and elastic, da e \/dt, pp and pp cross sections. 

Since the absorptive effects are included explicitly, the extracted values of the triple- Reggeon 
vertices are now much closer to the bare triple-Regge couplings. The value obtained, 



is about a factor of 3 larger than that found in the original analyses, which did not include 
absorptive corrections. Here, is the Pomeron-proton coupling. The new value of the cou- 
pling, gsp, is consistent with a reasonable extrapolation of the perturbative BFKL Pomeron 
vertex to the low scale region [22]. Note also that the resulting values of the 'PPP' and l PPK 
vertices allow a good description [23] of the HERA data [24] on inelastic J/ifi photoproduction, 
7P — > J/ip + Y, where the screening corrections are rather small. 

2.4 Multi-component s- and ^-channel model of soft processes 

Since the effects due to the triple-Pomeron vertex (8) are rather large, we must include all the 
multi- Pomeron diagrams, some of which are shown in Fig. 2(c). Why do we claim a large effect 
when A is only 0.2? The reason is that the contribution caused by such vertices is enhanced by 
the logarithmically large phase space available in rapidity. In particular, the total cross section 
of high-mass dissociation is roughly 7 of the form 



where A reflects the suppression of high-mass dissociation in comparison with elastic scattering 
and the Ins factor comes from the integration / dM 2 /M 2 ~ Ins. Thus actually we deal with 
the parameter Alns > 1 at collider energies. For each fixed rapidity interval the probability of 
high-mass dissociation (or, in other words, the contribution due to the triple-Pomeron vertex) 
is relatively small. However the cumulative effect in the complete interaction amplitude is 
enhanced by the large phase space available in rapidity. 

As a consequence, the contribution of the corresponding, so-called 'enhanced', diagrams, 
with a few vertices, is not negligible. Moreover, we expect that more complicated multi- Pomeron 
interactions, driven by the vertices (which describe the transition of n to m Pomerons), will 
affect the final result. Certainly[25], it is more reasonable to assume that oc \ n+m than 
to assume that g^ = for any n + m > 3. Thus, we need a model which accounts for the 
possibility of multi-Pomeron interactions (with arbitrary n and m). 

Here we follow the partonic approach of Refs. [26, 27]. While the eikonal formalism describes 
the rescattering of the incoming fast particles, the enhanced multi-Pomeron diagrams represent 

7 Hcre, for simplicity, we assume an essentially flat energy dependence, a ~ s e with elns < 1. The final 
equality in (9) can be deduced by taking the ratio of the couplings indicated in the Regge expressions in Fig. 1 . 



93P = ^9n, 



with A ~ 0.2, 



(8) 




(9) 
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Figure 3: (a) The ladder structure of the triple-Pomeron amplitude between diffractive eigen- 
states \4>i)i \4>k) of the proton; the rapidity y spans an interval to Y = Ins. (b) A multi- 
Pomeron diagram, (c) The "nested" ladder structure of the g\ multi-Pomeron vertex. 



the rescattering of the intermediate partons in the ladder (Feynman diagram) which describes 
the Pomeron-exchange amplitude. We refer to Fig. 3. The multi-Pomeron effects are included 
by the following equation describing the evolution in rapidity y of the opacity flk starting from 
the 'target' diffractive eigenstate \4>k)'- 

dn k (y,b) = e - AfW ,b)/ 2 e -An fc (»,6)/2 (aW^-Wm). (10) 



dy 

where y' = Ins — y. Let us explain the meanings of the three factors on the right-hand-side 
of (10). If only the last factor, (...)f4, is present then the evolution generates the ladder- 
type structure of the bare Pomeron exchange amplitude, where the Pomeron trajectory ap = 
1 + A + a't. The inclusion of the preceding factor allows for rescatterings of an intermediate 
parton c with the "target" proton k; Fig. 3(a) shows the simplest (single) rescattering which 
generates the triple-Pomeron diagram. Finally, the first factor allows for rescatterings with 
the beam i. In this way the absorptive effects generated by all multi-Pomeron diagrams are 
included, like the one shown in Fig. 3(b). There is an analogous equation for the evolution in 
rapidity y' of fli(y',b) starting from the 'beam' diffractive eigenstate \4>i)- The two equations 
may be solved iteratively. 

As we are dealing with elastic amplitudes we use e^ xn ^ 2 and not e~ xn . The coefficient A in 
the exponents arises since parton c will have a different absorption cross section from that of 
the diffractive eigenstates. Naively, we may assume that the states i, k contains a number 1/A 
of partons. The factors e~ xn ^ 2 generate multi-Pomeron vertices of the form 

g n m = nm X n+m - 2 g N /2 for n + m > 3 , (11) 

where a factor 1/n!, which comes from the expansion of the exponent, accounts for the identity 
of the Pomerons. The factors n{m) allow for the n(m) possibilities to select the Pomeron Q^(flfc) 
which enters the evolution (10) from the n(m) identical Pomerons. 

Even though A ~ 0.2 — 0.25, the role of factors e~ An / 2 is not negligible, since the suppression 
effect is accumulated throughout the evolution. For instance, if A <C 1 the full absorptive 
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correction is given by the product XQY/2, where the small value of A is compensated by the 
large rapidity interval Y . 

Strictly speaking nothing is known, either experimentally or theoretically, about the be- 
haviour of multi-Pomeron vertices, g 7 ^, at low scales. At large scales, for a large number of 
colours, N c , the leading contribution is given by diagrams like Fig. 3(a,c), which contain the 
g\ and g\ multi-Pomeron vertices respectively, where the Pomeron ladders are 'nested' within 
each other. Here the factor 1/n! arises from the time ordering (that is, from the ordering of the 
longitudinal momentum fractions) in the cells of each Pomeron adjacent to the parton c line; 
the lifetime of an inner Pomeron should be smaller than that for an outer Pomeron. Thus, in 
the large N c limit, the form (11) of the multi-Pomeron vertices, with A oc N c a s /ir, is motivated 
by perturbative QCD. At moderate values of N c , unfortunately there are other, more compli- 
cated, contributions. In particular, the 2 — » 2 Pomeron vertex g\ has a term proportional to 
a s /(Nc — l) 2 , that is to the first power of a s [28], and not the third power of a s as in (11). 
(Note that, in perturbative QCD, g^ oc a s as well.) However, this term is strongly suppressed 
by the colour factor 1/(N% — l) 2 = 1/64 and, most probably, such terms are not too important 
at low scales, when the coupling a s is not very small. 

So far, we have allowed multi-components in the s-channel via a multichannel eikonal. 
However, a novel feature of the model of Ref. [27] is that four different t-channel states are 
included. One for the secondary Reggeon (R) trajectory, and three Pomeron states (Pi, P 2 , P3) 
to mimic the BFKL diffusion in the logarithm of parton transverse momentum, ln(k t ) [29]. 
Recall that the BFKL Pomeron [30] is not a pole in the complex j-plane, but a branch cut. 
Here the cut is approximated by three t-channel states of a different size. The typical values of 
kt are kn ~ 0.5 GeV, k t 2 ~ 1.5 GeV and kt3 ~ 5 GeV for the large-, intermediate- and small-size 
components of the Pomeron, respectively. Thus (10) is rewritten as a four-dimensional matrix 
equation for Q% in t-channel space (a = Pi, P 2 , P3, R), as well as being a three-channel eikonal 
in diffractive eigenstate \(f>k) space. The transition terms, added to the equations, which couple 
the different t-channel components, are fixed by the properties of the BFKL equation. So, in 
principle, we have the possibility to explore the matching of the soft Pomeron (approximated by 
the large-size component Pi) to the QCD Pomeron (approximated by the small-size component 
P3). The key parameters which drive the evolution in rapidity are the intercepts 1 + A a and 
the slopes a' a of the t-channel exchanges. 

Section 3 will be devoted to the exclusive central production of different heavy systems 
A. That is, to processes of the type pp — > p + A + p. In Section 3.3, we will see that the 
multi-component Pomeron will allow us to obtain an estimate of the survival probability, S^ nh , 
of the rapidity gaps (denoted by the + signs) to so-called enhanced soft rescattering; that is, 
to rescattering on intermediate partons with different k t . Clearly enhanced rescattering will 
violate soft-hard factorisation of the process. Note that S^h depends on k t , which is driven by 
the scale of the central hard subprocess. We need to introduce components of the Pomeron of 
different size to be able to calculate S en h for different central exclusive processes, such as A = 
Higgs, dijet, 77 or Xc- These enhanced rescattering effects are in addition to the survival of the 
gaps to eikonal soft rescattering, which preserves soft-hard factorisation. 
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-t (GeV 2 ) 

Figure 4: The t dependence of the elastic pp cross section [27]. The dashed and dotted lines 
are the contributions from the elastic scattering of the largest size (i — 1) and the smallest size 
(i = 3) diffractive components. 

2.5 Description of the 'soft' data and predictions for the LHC 

The number of parameters in the model of Ref. [27] is too large to perform a straightforward 
X 2 fit of the data. Instead, the majority of the parameters are fixed at physically reasonable 
values, and it was demonstrated that all the features of the available data on diffractive cross 
sections, cr tot , da c \/dt, <j]^ m , da S T>/dM 2 , in the CERN-ISR to Tevatron range are reproduced. 

Fig. 4 shows the quality of the description [27] of the data for the elastic differential cross 
section. We also present in Fig. 4 the prediction for differential elastic cross section at the LHC 
energy y/s = 14 TeV. Recall that [27] uses a three-channel eikonal. That is i, k — 1, 2, 3. It is 
interesting to note that the contribution to the cross section arising from the scattering of the 
two large-size diffractive eigenstates, (i — 1) x (i — 1), already has a diffractive dip at —t = 0.2 
GeV 2 . However, after the contributions from all possible combinations i x k are summed up, 
the prediction has no dip up to —t = 0.5 GeV 2 . 
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To describe the high-energy behavior of the total cross section, A a = a a (0) — 1 = 0.3 is taken 
for each of the three components of the Pomeron. These Pomeron intercepts are consistent with 
resummed NLL BFKL, which gives u ~ 0.3 practically independent of the scale k t [31]. The 
slopes of the Pomeron trajectories are driven by the transverse momentum associated with 
the particular component a; in fact, a' oc l/fcf. The data require a' Pi = 0.05 GeV~ 2 for the 
large-size Pomeron component, so in the model we take a'p 2 = 0.05/9 GeV~ 2 for the second 
component and a' P = for the smallest-size component. For the secondary Reggeon trajectory 
the model takes a' R = 0.9 GeV~ 2 , and q;_r(0) = 0.6. The 'bare' intercept is a bit larger than 
|, since the final effective intercept is reduced by the absorptive corrections included in the 
evolution equation. The description of the total cross section data are shown in Fig. 5(a). 
The screening corrections arising from the 'enhanced' multi-Pomeron diagrams, that is from 
the high-mass dissociation, slow down the growth of the cross section with energy. Thus, the 
model predicts a relatively low total cross section at the LHC - a to t(LHC) ~ 90 mb 8 . 

The value of the parameter A, which controls the cross section of high-mass dissociation in 
the small £ (that is, large y) region, was found to be A ~ 0.25. The dependence of the cross 
section for high-mass dissociation, C,d 2 a/dtdC,, on £ = M 2 /s is compared with the Tevatron 
CDF data [34, 35] in Fig. 6. The results without the ixixP term are shown by the dashed lines. 
We also show in Fig. 6(a) (by the dotted line at small £) the prediction for the LHC energy. 

The energy behaviour of the cross sections are shown in Table 1 and Fig. 5. We do not 
quote the values of the cross section for double diffractive dissociation, <7 DD . The model of [27] 
will give results similar to those in Table 2 of Ref. [26] . These values of <t D d are in excellent 
agreement with the cross sections observed at the Tevatron. 



energy 


Ctot 


Cel 


_lowM 
°SD 


highM 
°~SD 


tot 
°SD 


1.8 


73.7 


16.4 


4.1 


9.7 


13.8 


14 


91.7 


21.5 


4.9 


14.1 


19.0 


100 


108.0 


26.2 


5.6 


18.6 


24.2 



Table 1: Cross sections (in mb) versus collider energy (in TeV) [27]. 

The values of Og^ quoted in Table 1 look, at first sight, too large, when compared with 
the value 9.46 ± 0.44 mb given by CDF [34]. However the CDF value does not include the 
secondary Reggeon (RRP) contribution, denoted as a 'non- diffractive' component of 2.6 ± 0.4 
mb. Moreover, the trigger used to select the diffractive dissociation events rejects part of the 
low-mass proton excitations. Taking these absences into account, there is no contradiction 
between the model prediction and the CDF data. Furthermore, note that in the region where 

8 Low values, ~90 mb, are also predicted by other models of 'soft' interactions which include absorptive 
effects [32, 33]. 
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Figure 5: The energy dependence of the total (a), elastic and diffractive dissociation (b) pp 
cross sections and the cross sections of dissociation to a fixed M 2 = £s state (c); (d) the parton 
multiplicity (solid lines) and the number of 'colour tubes' (dashed) produced by the Pomeron 
components of different size. The figure is taken from [27]. 



the CDF detector efficiency and resolution are good, our model gives a good description of the 
differential cross sections which were actually measured, see Fig. 6. 

As we have a detailed model for high energy soft processes, it would appear to be possible 
to predict the multiplicity distribution at the LHC However, although some general features 
can be predicted, unfortunately, we cannot reliably calculate the multiplicity distributions of 
secondary hadrons, since a non-negligible fraction of the final hadrons may be produced via 
the fragmentation of minijets. We can only discuss the multiplicity distribution at the partonic 
level. 

The mean number of the (t-channel) ladders of the type a produced in the collision of % 
and k Good- Walker eigenstates can be calculated. This quantity may be considered as the 
mean number of colour tubes of type a produced in the proton-proton interaction. To obtain 
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Figure 6: The model description [27] of the data for the cross section for high-mass dissociation 
versus £ for -t = 0.05 GeV 2 at ^/s = 1800 GeV and 546 GeV [34, 35]. The dashed lines are 
the predictions without the unP contribution. The dotted curve at small £ is the prediction 
for the LHC. 



the number of partons created by the ladder 'a' at rapidity y, we have to allow for the parton 
density p a {y)- The results are shown in Fig. 5(d). The main growth in multiplicity, as we go 
from Tevatron to LHC energies, is due to the small size ('QCD') Pomeron component, which 
produces particles with typically pt ~ 5 GeV. There is essentially no growth in multiplicity at 
small p t . This simply confirms the trend that has been observed through the CERN-ISR to 
Tevatron energy range, see the data points in Fig. 7. 

In other words, starting with the same intercepts (A = 0.3) the contribution of the large- 
size component after the absorptive correction becomes practically flat, while the small-size 
contribution, which is much less affected by the absorption, continues to grow with energy. As 
mentioned above, such a behaviour is consistent with the experiment (see Fig. 7) where the 
density of low k t secondaries is practically saturated while probability to produce a hadron with 
a large (say, more than 5 GeV) transverse momentum grows with the initial energy. To obtain 
a qualitative feel for the expected behaviour at the hadronic level, we show the estimates of 
[27] at the Tevatron and LHC energies in Fig. 7, where the horizontal lines indicate the typical 
p t interval associated with each Pomeron component. 



3 Exclusive processes at hadron colliders 

In this Section we discuss exclusive processes of the type pp — > p + A + p where the + signs 
denote large rapidity gaps. The mechanism for the process in sketched in Fig. 8. Although, 
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Figure 7: The plot is from Ref. [36]. The horizontal lines, which are superimposed, are the 
predictions of [27] at the Tevatron and LHC energies; the three p t ranges correspond to the 
large-, intermediate- and small-size components of the Pomeron. 



below, we will consider several different heavy systems A, the main motivation is the possibility 
to use the process pp — > p + H + p to detect and study one or more Higgs bosons at the LHC. 

3.1 Advantages of the exclusive Higgs signal 

We consider the mass range Mh ^ 140 GeV, where the dominant decay mode is H — > bb. The 
exclusive process has a small cross section; cr exc i ~ 10~ 4 0"*°^ at the LHC energy. Nevertheless, 
it has the following advantages: 

• The mass of the Higgs boson (and in some cases the width) can be measured with high 
accuracy (with mass resolution a(M) ~ 1 GeV) by measuring the missing mass to the 
forward outgoing protons, provided that they can be accurately tagged some 400 m from 
the interaction point [12]. 
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• The leading order bb QCD background is suppressed [37] by the P-even J z = selection 
rule 9 , where the z axis is along the direction of the proton beam. Indeed, at LO, this 
background vanishes in the limit of massless b quarks and forward outgoing protons. 
Therefore, one can observe the Higgs boson via the main decay mode H — > bb, which 
otherwise is very challenging to measure in inclusive production due to the overwhelming 
QCD background. Moreover, a measurement of the mass of the decay products must 
match the 'missing mass' measurement. 

• The quantum numbers of the central object (in particular, the C- and P-parities) can be 
analysed by studying the azimuthal angle distribution of the tagged protons. Due to the 
selection rules, the production of ++ states is strongly favoured. 

• There is a clean environment for the exclusive process - the soft background is strongly 
suppressed [38]. 

• Extending the study to BSM Higgs bosons [39] - [45]. For example, there are regions of 
MSSM parameter space were the signal is enhanced by a factor of 10 or more, while the 
background remains unaltered [39, 41, 43]. Moreover, there are domains of parameter 
space where MSSM Higgs boson production via the conventional inclusive processes is 
suppressed whereas the exclusive signal is enhanced, and even such, that both the h and 
H ++ bosons may be detected [41]. 

The main background to the exclusive H — > bb signal comes from (i) irreducible QCD 
gg pp — > bb events, (ii) gluons mimicking b jets and (iii) the \ J Z \ =2 contribution arising from 
non-forward going protons. The superscript PP is to indicate that the active gluons originate 
within overall colourless hard Pomeron exchanges. Recently [46], the NLO contribution to (i) 
has been calculated and reduces the irreducible background by a factor of 2 or more. Radiation 
from the screening gluon in Fig. 8 is numerically small [47]. In summary, it should be possible 
to achieve a H — > bb signal-to-background ratio of 0(1), and much greater in some BSM Higgs 
scenarios [41, 42, 45]. 

The downside is that the exclusive cross section is small: about 3 fb for a 120 GeV SM 
Higgs at the LHC. Allowing for acceptance cuts, b jet recognition etc., it corresponds to only 
a handful, or so, observable events for an integrated luminosity of 30 fb -1 . Moreover, there 
are still some experimental issues to be settled; in particular, very good timing is necessary to 

9 J z — originates from s-channel helicity conservation of the forward protons if their transverse momenta 
Pa = 0, with i = 1,2. The admixture of the \ J Z \ = 2 amplitudes is governed by the product (p t i • r) (p t 2 • 
r) ~ (pf) /Qt, where the transverse size of the digluon Pomeron r ~ 1 /Qt, and Qt is the transverse momentum 
in the gluon loop in Fig. 8. Thus, the \J Z \ = 2 contribution to the cross section is suppressed by a factor p\jQ\- 
The presence of the Sudakov form factor provides the infrared stability of the Qt integral over the gluon loop. 
Typically, the main contribution comes from the saddle-point region of the loop integral, namely Q\ ~ 4 GeV 2 
for the exclusive production of a scalar Higgs of mass Mh = 120 GeV at the LHC energy, and so ensures that 
J z = is a very good approximation. The P even selection reflects the symmetry between the active gluons 
emitted from the two incoming protons. 
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Figure 8: The mechanism for the exclusive process pp — ► p + A + p, with the eikonal and 
enhanced survival factors shown symbolically. The thick lines on the Pomeron ladders, either 
side of the subprocess (gg — > A), indicate the rapidity interval Ay where enhanced absorption 
is not permitted, see Section 3.3. 



identify an exclusive event from among the multiple interactions per bunch crossing which will 
occur at high LHC luminosities. 



3.2 Calculation of the exclusive Higgs signal 

The calculation of the exclusive production of a heavy system is an interesting mixture of soft 
and hard QCD effects. Here we concentrate on exclusive scalar Higgs production, pp — > p+H+p. 
The basic mechanism is shown in Fig. 8. The t-integrated cross section is of the form 



S 2 

a — — 

5 2 



(12) 



where B/2 is the t-slope of the proton- Pomeron vertex, and N is given in terms of the H — > gg 
decay width. The probability amplitudes, f g , to find the appropriate pairs of i-channel gluons 
(xi,x[) and (x2,x' 2 ), are given by the skewed unintegrated gluon densities at a hard scale 
ji ~ Mff/2. Since (x' ~ Qt/Vs) "C (x ~ Mh/s/s) 1, it is possible to express f g (x,x' ,Q 2 , fi 2 ), 
to single log accuracy 10 , in terms of the conventional integrated density g(x), together with a 
known Sudakov suppression factor T, which ensures that the active gluons do not radiate in 
the evolution from Q t up to the hard scale \x ~ Mjj/2, and so preserves the rapidity gaps. The 
factor T ensures that the integral is infrared stable, and may be calculated by perturbative 
QCD. 

The factor S 2 in (12) is the probability that the secondaries, which are produced by soft 
rescattering do not populate the rapidity gaps. It is the price we have to pay for an exclusive 



10 This can be achieved, to a good approximation, by setting the lower limit of the k t integration in the 
Sudakov form factor equal to Qt [48], based on the BFKL equation; and the upper limit fj, = 0.62Mh- The 
factor 0.62 follows from an exact calculation of the one-loop contribution [39]. 
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process. As written, the cross section assumes soft-hard factorization. It assumes that the 
survival factor, denoted by in Fig. 8 and calculated from a model of soft interactions, does 
not depend on the structure of the perturbative QCD amplitude embraced by the modulus 
signs in (12). Actually the situation is more complicated, see, for instance, [16, 48]. 



3.3 Rapidity gap survival 



The gap survival factor caused by eikonal rescattering of the diffractive eigenstates [19], for a 
fixed impact parameter b, is 



EH 2 \a k \ 2 M tk (b) exp(-Q^(s,b)/2) 

i,k 



EH 2 W\ 2 M ik (b) 

i,k 



(13) 



where ^^(s, b) is the total opacity of the ik interaction, and the a^'s occur in the decomposition 
of the proton wave function in terms of diffractive eigenstates \p) = J2i a «| </>«)■ The total opacity 
has the form Vt%(y)Q%(y') integrated over the impact parameters bi, b2 (keeping a fixed impact 
parameter separation b = bi — b2 between the incoming protons) and summed over the different 
Pomeron components a. Recall y' = Y — y = Ins — y, see Fig. 3. The exact shape of the matrix 
element M.,^ for the hard subprocess gg — > H in b space and the relative couplings to the 
various diffractive eigenstates i, k should be addressed further. 

One possibility is to say that the b dependence of M. should be, more or less, the same as 
for diffractive J/ip electroproduction (j+p — > J/ip+p), and the coupling to the \4>i) component 
of the proton should be proportional to the same factor % as in a soft interaction. This leads 
to 

M ik oc 7 i7fc exp(-b 2 /4B) (14) 

with t-slope B ~ 4 GeV~ 2 [49]. The resulting "first look" predictions obtained using the 'soft' 
model of [27], for the exclusive production of a scalar 120 GeV Higgs at the LHC, are shown 
in Fig. 9. After we integrate over b, we find that the survival probability of the rapidity gaps 
in pp — > p + H + p to eikonal rescattering is (S' 2 ik )=0.017, with the Higgs signal concentrated 
around impact parameter b = 0.8 fm. Expressing the survival factors in this manner is too 
simplistic and even sometimes misleading, for the reasons we shall explain below; nevertheless 
these numbers are frequently used as a reference point. 

As indicated in Fig. 8, besides eikonal screening, S^, caused by soft interactions between 
the protons, we must also consider so-called 'enhanced' rescattering, S^nh, which involves inter- 
mediate partons. Since we have to multiply the probabilities of absorption on each individual 
intermediate parton, the final effect is enhanced by the large multiplicity of intermediate par- 
tons. Unlike S^b), the enhanced survival factor S^ nh (h) cannot be considered simply as an 
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Figure 9: A "first look" at the impact parameter dependence of the signal for 120 GeV Higgs 
production at the LHC after including an eikonal rescattering correction. 

overall multiplicative factor. The probability of interaction with a given intermediate parton 
depends on its position in configuration space; that is, on its impact parameter b and its mo- 
mentum k t . This means that S en h simultaneously changes the distribution of the active partons 
which finally participate in the hard subprocess. It breaks the soft-hard factorization of (12). 

Do we anticipate that S en ^ will be important? Working at LO (of the collinear approxima- 
tion) we would expect that effect may be neglected. Due to strong /c r ordering the transverse 
momenta of all the intermediate partons are very large (i.e. the transverse size of the Pomeron 
is very small) and therefore the absorptive effects are negligible. Nevertheless, this may be not 
true at a very low x, say x ~ 10~ 6 , where the parton densities become close to saturation and 
the small value of the absorptive cross section is compensated by the large value of the parton 
density. Indeed, the contribution of the first enhanced diagram, which describes the absorption 
of an intermediate parton, was estimated in the framework of the perturbative QCD in Ref.[50]. 
It turns out that it could be quite large. On the other hand, such an effect does not reveal 
itself experimentally. The absorptive corrections due to enhanced screening must increase with 
energy This is not observed in the present data (see [51] for a more detailed discussion). One 
reason is that the gap survival factor S 2 ik already absorbs almost the whole contribution from 
the centre of the disk. The parton essentially only survives eikonal rescattering on the periph- 
ery; that is, at large impact parameters b. On the other hand, on the periphery, the parton 
density is rather small and the probability of enhanced absorption is not large. This fact can 
be seen in Ref. [52]. There, the momentum, Q s , below which we may approach saturation, was 
extracted from the HERA data in the framework of the dipole model. (Q s is the inverse size of 
the dipole for which absorptive corrections become sizeable.) Already at b = 0.6 fm the value 
of Q 2 S is such that Q 2 S < 0.3 GeV 2 for the relevant values of x, namely x ~ 10~ 6 . However, 
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in the perturbative QCD calculations of Refs. [10, 53, 54] the infrared cutoff Qt > Qo = 0.85 
GeV was introduced in order not to have uncontrollable uncertainties in the parton distribution 
functions coming from the non-perturbative domain. 

Now, the model of Ref. [27], with its multi-component Pomeron, allows us to calculate the 
survival probability of the rapidity gaps, to both eikonal and enhanced rescattering. Recall 
that the evolution equations in rapidity (like (10)) have a matrix form in aa' space, where 
a = 1,2,3 correspond to the large-, intermediate- and small-size components of the Pomeron. 
We start the evolution from the large component P±, and since the evolution equations allow 
for a transition from one component to another (corresponding to BFKL diffusion [29] in 
lnk t space), we determine how the enhanced absorption will affect the high-/^ distribution in 
the small-size component P3, which contains the active gluon involved in forming the Higgs. 
Moreover, at each step of the evolution the equations include absorptive factors of the form 
e -A(ng+n?)/2_ gy solving the equations with and without these suppression factors, we could 
quantify the effect of enhanced absorption. However, there are some subtle issues here, see 
[16]. First, since we no longer have soft-hard factorization, we must first specify exactly what 
is included in the bare hard amplitude. 

Another relevant observation is that the phenomenologically determined generalised gluon 
distributions are usually taken at pt = 0, and then the observed "total" cross section is cal- 
culated by integrating over p t of the recoil protons assuming an exponential behaviour e~ Bp * ; 



However, the total soft absorptive effect changes the p t distribution in comparison to that for 
the bare cross section determined from perturbative QCD. Moreover, the correct p t dependence 
of the matrix element M. of the hard gg — > H subprocess does not have an exponential form. 
Thus the additional factor introduced by the soft interactions is not just the gap survival S 2 , 
but rather S 2 {p 2 ) 2 , where the square arises since we have to integrate over the pt distributions of 
two outgoing protons. Indeed in all the previous calculations the soft prefactor had the form 11 
S 2 /B 2 . Note that, using the model of Ref. [27], we no longer have to assume an exponential b 
behav! iour of the matrix element. Now the b dependence of Ai(h) is driven by the opacities, 
and so is known. Thus we present the final result in the form S 2 (p 2 ) 2 . That is, we replace 
S 2 /B 2 in (12) by S 2 {p 2 } 2 . So if we wish to compare the improved treatment with previous 
predictions obtained assuming B = 4 GeV -2 we need to introduce the "renormalisation" factor 
((p 2 )B) 2 . The resulting (effective) value is denoted by S 2 S . 

11 At larger impact parameter b the absorption is weaker. Hence the value of S 2 increases with the slope B. 
It was shown that the ratio S 2 /B 2 is approximately stable for reasonable variations of B [53]. 



that is 




(15) 



where 




(16) 
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Before we do this, there is yet another effect that we must include. We have to allow for a 
threshold in rapidity [55, 51, 56]. The evolution equation for (10), and the analogous one 
for Qf, are written in the leading ln(l/x) approximation, without any rapidity threshold. The 
emitted parton, and correspondingly the next rescattering, is allowed to occur just after the 
previous step. On the other hand, it is known that a pure kinematical i m i n effect suppresses 
the probability to produce two partons close to each other. Moreover, this t min effect becomes 
especially important near the production vertex of the heavy object. It is, therefore, reasonable 
to introduce some threshold rapidity gap, Ay, and to compute S 2 nh only allowing for absorption 
outside this threshold interval, as indicated in Fig. 8. For exclusive Higgs boson production at 
the LHC 12 , the model gives S 2 ff = 0.004, 0.009 and 0.015 for Ay = 0, 1.5 and 2.3 respectively 
[16]. For Ay = 2.3 all the NLL BFKL corrections [31] may be reproduced by the threshold 
effect [55, 51, 56, 58]. 

Furthermore, Ref. [16] presents arguments that 

(S 2 eS ) = 0.015 + _ f 05 (17) 

should be regarded as a conservative (lower) limit for the gap survival probability in the ex- 
clusive production of a SM Higgs boson of mass 120 GeV at the LHC energy of y/s = 14 TeV. 
Recall that this effective value should be compared with S 2 obtained using the exponential 
slope B = 4 GeV -2 . The resulting value for the cross section is, conservatively, 

a(pp^p + H + p) ~ 2-3 fb, (18) 

with an uncertainty 13 of a factor of 3 up or down, see also [39, 48]. In Section 3.6 we explain 
why the values (17) and (18) should be regarded as conservative. 

12 A very small value S^ nh = 0.063 is claimed in [33], which would translate into an extremely small value of 
S^ s = 0.0235 x 0.063 = 0.0015. There are many reasons why this estimate is invalid. In this model the two- 
particle irreducible amplitude depends on the impact parameter b only through the form factors of the incoming 
protons. The enhanced absorptive effects (which result from the sum of the enhanced diagrams) are the same 
at any value of b. Therefore, the enhanced screening effect does not depend on the initial parton density at 
a particular impact parameter 6, and does not account for the fact that at the periphery of the proton, from 
where the main contribution comes (after the S e ik suppression), the parton density is much smaller than that in 
the centre. For this reason the claimed value of S* 2 nh is much too small. Besides this lack of k t <-> b correlation, 
the model has no diagrams with odd powers of g$p. For example, the lowest triple-Pomeron diagram is missing. 
That is, the approach does not contain the first, and most important at the lower energies, absorptive correction. 
Next, recall that in a theory which contains the triple-Pomeron coupling only, without the four-Pomeron term 
(and/or more complicated multi-Pomcron vertices), the total cross section decreases at high energies [57]. On 
the other hand, the approximation used in [33] leads to saturation (that is, to a constant cross section) at very 
high energies. In other words, the approach is not valid at high energies. This means that such an approximation 
can only be justified in a limited energy interval; at very high energies it is inconsistent with asymptotics, while 
at relatively low energies the first term, proportional to the first power of the triple-Pomeron coupling g^p, is 
missing. Finally, the predictions of the model of [33] have not been compared to the CDF exclusive data of 
Section 3.4. 

13 Bcsides the uncertainty arising from that on S^ s , the other main contribution to the error comes from that 
on the unintegrated gluon distributions, f g , which enter to the fourth power [48, 43]. 
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Figure 10: The mechanism for the exclusive processes observed by CDF at the Tevatron. The 
survival probabilities of the rapidity gaps are not shown in the sketches. 

3.4 Exclusive processes observed at the Tevatron 

Exclusive diffractive processes of the type pp — > p + A+p have already been observed by CDF 
at the Tevatron, where A = 77 [59] or dijet [60] or Xc [61]. As the sketches in Fig. 10 show, 
these processes are driven by the same mechanism as that for exclusive Higgs production, but 
have much larger cross sections. They therefore serve as "standard candles" . 

CDF observe three candidate events for pp — > p + 77 + p with Ej, > 5 GeV and |^ 7 | < 1 
[59]. Two events clearly favour the 77 hypothesis and the third is likely to be of 7r°7r° origin. 
The predicted number of events for these experimental cuts is O.SIqJ [62], giving support to the 
theoretical approach used for the calculation of the cross sections for exclusive processes. In 
fact, central exclusive two- and three-jet production offers more detailed tests of the theoretical 
formalism, see [48]. 

Especially important are the recent CDF data [60] on exclusive production of a pair of high 
E T jets, pp — > p+jj+p. As discussed in [63, 9, 10] such measurements could provide an effective 
gg pp 'luminosity monitor' just in the kinematical region appropriate for Higgs production. 
The corresponding cross section was evaluated to be about 10 4 times larger than that for the 
production of a SM Higgs boson. Since the exclusive dijet cross section is rather large, this 
process appears to be an ideal 'standard candle'. A comparison of the data with analytical 
predictions [63, 10] is given in Fig. 11. It shows the E^ m dependence for the dijet events 
with Rjj = Mdijet/Mpp > 0.8, where Mpp is the invariant energy of the incoming Pomeron- 
Pomeron system. The agreement with the theoretical expectations [9, 10] lends credence to the 
predictions for the exclusive Higgs production. 

In particular, these CDF dijet data clearly demonstrate the suppression due to the Sudakov 
factor T which occurs in the evaluation of the unintegrated gluon distribution, f g , needed to 
predict the exclusive cross section, (12). On dimensional grounds we would expect da jdE\ oc 
\jE\. This behaviour is modified by the anomalous dimension of the gluon and by a stronger 
Sudakov suppression with increasing Ex- Already the existing CDF exclusive dijet data [60] 
exclude predictions which omit the Sudakov effect. In Section 3.7 we describe how observation 
of central exclusive three-jet, as well as two-jet, production at the LHC offers more detailed 
tests of f g . 
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Figure 11: The cross section for 'exclusive' dijet production at the Tevatron as a function E™ m 
as measured by CDF [60]. The data integrated over the domain Rjj = M&^/Mpp > 0.8 
and Et > E™ m - A jet cone of R < 0.7 is used. The curves are the pure exclusive cross 
section calculated [63, 10] using the CDF event selection. The solid curve is obtained [64] by 
rescaling the parton (gluon) transverse momentum pt to the measured jet transverse energy Et 
by Et = 0.8pr- The dashed curve assumes Et = 0.75pT- The rescaling procedure effectively 
accounts for the hadronization and radiative effects, and for jet energy losses outside the selected 
jet cone. 




CDF [61] have also recently observed exclusive Xc production, pp — > p + Xc+P, via the decay 
chain Xc — > J /4> + 7 — ► P>~1- They collected 65 ± 10 signal events, with limited M(J /iffj) 
resolution. Assuming the dominance of Xc(0) production, this corresponds to 

da(xc)/dy\ y=0 = 76 ± 14 nb. (19) 

Strictly speaking, due to the low scale associated with exclusive Xc production, only order- 
of-magnitude estimates can be made for the cross section. Ref. [53] assumes a perturbative 
contribution coming from integrating round the gluon loop in Fig. 8 for Q t > 0.85 GeV, with 
the remaining imfrared contribution estimated non-perturbatively. This gave a prediction for 
da(xc(0 ++ ) /dy\ y= o of 90 nb, 14 with | coming from the non-perturbative region and | from the 
perturbative domain. The rapidity gap survival factor was taken to be (S 2 ) = 0.07, assuming 
only eikonal rescattering. However, for the relatively light Xc(0), the available rapidity interval 
for the enhanced suppression is large. In this case the inclusion of enhanced rescattering may 

14 The value 130 nb obtained in [53] has been changed to 90 nb to take into account the revised value of the 
total Xc(0) width of 10.2 MeV as compared to the value in the RPP(2002) which was 1.45 times higher. 
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cause a large difference - reducing the exclusive cross section by up to |. The resulting cross 
section of 30 nb, which includes S^, is below the value given in (19). 

However, in the conditions of the experiment [61], where the transverse momenta of the 
forward outgoing protons are not measured, the situation for Xc m &y be more complicated, 
due to possible contributions from the higher spin % c (l ++ ) and x c (2 ++ ) states. Recall that 65 
events with two large rapidity gaps were observed via the Xc J / ^7 decay mode, and that the 
Xc(0) branching fraction to this mode is very low. Explicitly, the Xc(0), Xc(1)j Xc(2) branching 
fractions to J/ip'y are 0.013, 0.36 and 0.20 respectively. Also, the bare cross sections for Xc(l) 
and Xc(2) production are suppressed relative to x c (0) production by factors of {p^)/M 2 and 
(pi) jQ\ respectively, where p t is the transverse momentum of the outgoing protons and Q t 
is that round the gluon loop. Numerically this leads to a suppression by about a factor of 
30 — 100, depending on the choices of the global partons and of the infrared cutoff that are 
used in the computation. Part of this suppression is compensated by a larger gap survival 
probability for x c (l) an d Xx(2), since, due to their spin structure, they are produced more 
peripherally. Thus, finally, we estimate a production ratio x c (0)/x c (l) ~ 10 — 40, with more 
or less the same suppression ratio found for x c (2) 15 . In summary, since the energy resolution 
of the CDF observation is insufficient to separate the three Xc states, then, allowing for the 
J/ipj branching fractions, we may expect more or less comparable contributions from each of 
the states 16 . Therefore even the conservative expectation of 30 nb for Xc(0) is not inconsistent 
with (19). 

To make further progress, it would be instructive to observe central exclusive Xc production 
via the or KK decay channels, see [53]. Recall that the rnr or KK decay modes of Xc(0) 
have a branching fraction of about 1%, while these decay channels are forbidden for XciX), and 
suppressed by about a factor of 5 for Xc(2) relative to x c (0); rather than enhanced relative 
to Xc(0), as was the case for the J/ipl channel. Alternatively we could study the angular 
distribution of the J/ipj system in order to distinguish the spins of the parent Xc- 

As a final general comment, we note that it is sometimes stated that enhanced diagrams, 
which break soft-hard factorization, generate an extremely small gap survival factor, S^nh [33, 
66]. However, the analysis of [16] shows the following hierarchy of the size of the gap survival 
factor to enhanced rescattering 

SZ C (M H > 120 GeV) > SjT ron (ir, E T > 5 GeV) > SjT ron (Xc), (20) 

which reflects the size of the various rapidity gaps of the different exclusive processes. The fact 
that 77 and Xc events have been observed at the Tevatron confirms that there is no danger that 
enhanced absorption will strongly reduce the exclusive SM Higgs signal at the LHC energy. 

15 Recall that these predictions at such low scales have a large uncertainty. Moreover, the NLO corrections 
could be quite sizeable, and require further detailed studies. Note, also, that the relative number of events 
where the forward protons dissociate is larger for Xc(l) and Xc(2) than for Xc(0). 

16 The importance of x c (l) production has been recently emphasized in Ref. [65]. 
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3.5 Dependence of the survival factors on collider energy 

It is relevant to ask how the survival factors, S 2 ik and S 2 nh , change in size in going from the 
Tevatron to the LHC energy. The energy dependence of the eikonal gap survival S 2 ik is driven 
by the ratio a to t/B e \, which controls the behaviour of the opacity Q. Here, B e \ is the slope of 
the elastic differential cross section. As we have seen, the growth of the total cross section is 
not expected to be large, and is partly compensated by the increase in the value of the slope 
B e \. Indeed, for exact "geometric scaling", for which a tot oc R 2 (where R is the interaction 
radius), the function Q(b) does not depend on energy at all. Thus we expect a rather weak 
(only logarithmic) decrease of the eikonal survival probability with energy. 

The enhanced survival factor S 2 nh may decrease with energy due to 

(i) the increased size of the rapidity interval (oc ln(s/M 2 )) available for screening corrections, 

(ii) the increase in the parton density as x ~ 1/s decreases, 

where we are concerned with the part of the amplitude denoted by S cnh in Fig. 8. We discuss 
these possible s dependences in turn. 

At fixed M 2 , the available rapidity interval grows with s and this should lead to a decrease 

of 

$enh- Recall that M is the mass of the centrally produced system. However, no sizeable 
variation with energy is observed in present data, say, in leading neutron production at HERA. 
This phenomenological fact indicates that the role of S 2 nh factor should be relatively small (see 
[51] for a detailed discussion). 

The values of x ~ relevant at LHC energies are quite small. At leading order (LO) 
the gluon density increases with 1/x. So, at first sight, it appears that this may lead to the 
Black Disk Regime where the enhanced absorptive corrections will cause the cross section for 
exclusive production to practically vanish. However, NLO global parton analyses show that at 
relatively low scales (k 2 = /i 2 ~ 2 — 4 GeV 2 ) the gluons are flat for x < 1CT 3 — 10~ 4 (or even 
decrease when x — * 0). Recall that the contribution of enhanced diagrams from larger scales 
decrease as 1/k 2 (see [50]). The anomalous dimension 7 < 1/2 is not large and so the growth of 
the gluon density, xg(x, fx 2 ) oc (/i 2 ) 7 , cannot compensate the factor 1/A; 2 . Therefore the whole 
enhanced diagram contribution should be evaluated at rather low scales where the NLO gluon 
is flat 17 in x. Thus there is no reason to expect a strong energy dependence of S 2 nh . 

3.6 Why the estimates are conservative 

It is useful to list the reasons why we regard the estimates of the rapidity gap survival factor 
(17) and the exclusive cross section (18) as conservative lower limits. 

17 Note that the flat x-behaviour of NLO gluons allow the justification of the inequalities in (20). 
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First, we should not use the whole irreducible amplitude Q^b) when calculating the ab- 
sorption (as, for example, in exp(— r2^ t (b)/2) of (13)), but only that part of it which corresponds 
to having some rescattering which populates the gap in the chosen rapidity interval. The com- 
putations of [27, 16] neglect this effect. This means that the gap survival probabilities, and the 
true cross sections of diffractive dissocation, should be a bit larger than the predictions given 
earlier. 

Second, a detailed study [16] of the b dependence of S^ nh , obtained from the soft model of 
[27], is instructive. Assuming that all the partons are distributed homogeneously, we would 
expect S^ nh — > 1 at the (large) values of b occurring in the periphery of the pp interaction. 
However, the results show only a weak tendency to increase at large b. This may be explained 
by the fact that in the 'soft' model of [27] (with its small value of a' P ), S^ nh comes mainly from 
u hot-spots" in which many individual intermediate partons are concentrated within small b 
domains. In other words, most of enhanced absorption occurs within the same parton shower, 
and is due to secondaries produced during the evolution with practically the same impact 
parameter b. On the other hand, a detailed analysis [52] of HERA data shows that the value 
of the saturation scale Q s (b) decreases rapidly with increasing b. (Recall that Q s is the inverse 
size of the dipole for which absorptive corrections become important.) This indicates that the 
parton-parton correlations are too strong in the present model. Therefore, we consider our 
results are close to the maximum possible gap suppression. 

Third, the partonic nature of the diffractive eigenstates is not specified in the soft model of 
[27]. When BFKL is introduced, it is implied that the partons are gluons, and so the screening 
exponents were chosen assuming gluons in both Pomerons. However, at NLO the gluon density 
in the relevant domain is approximately flat in x. In terms of QCD the increase of the parton 
density with decreasing x, at the relevant low scales, may reflect the growth of the NLO sea 
quark density with decreasing x. On the other hand, the central exclusive amplitudes of interest 
are driven by gluon-gluon fusion. If there is screening it should be caused by the (NLO) sea- 
quark contribution in these amplitudes. Then the effective colour factor is smaller than in the 
model of [27] where BFKL (that is, pure gluonic colour coefficients) is used to evaluate the 
absorptive effects. Hence, again, we expect a larger value of S^. 

Finally, it is worth mentioning that in deriving predictions for the cross sections of central 
exclusive processes, such as (18), rather conservative choices of parton distributions were used 
[10, 39, 53]. Larger predictions will result if the gluon densities are found to be even a little 
larger, recall the (xg) A dependence. 

We note from Section 3.4 that, on average, the predictions lie a little below the present 
exclusive data obtained by CDF at the Tevatron. However this is in line with the above 
comments. It will be informative to see how the data/theory comparisons hold up with the 
improved experimental statistics which will come from the Tevatron and, later on, from the 
LHC. 
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Figure 12: (a) W production with 2 gaps, (b) Inclusive W production, (c) Z production with 
2 gaps. 



3.7 Early LHC probes of exclusive production 

Recall that the uncertainties associated with predictions for an exclusive process are potentially 
not small. Each stage has its own uncertainties. Therefore, it is important to perform checks 
of the approach using processes with appreciable cross sections that will be experimentally 
accessible in the first data runs of the LHC; that is with integrated luminosities in the range 
100 pb _1 to 1 fb _1 . In fact, it is possible identify processes where the different ingredients of 
the formalism used to calculate central exclusive production can be tested experimentally, more 
or less independently [48]. 

To probe the gap survival factor to eikonal rescattering, 5f ik , we need a process with a bare 
cross section that can be calculated reliably. Good candidates are the production of W or Z 
bosons with rapidity gaps on either side [48]. In the case of W+gaps' production the main 
contribution comes from the diagram of Fig. 12(a). One gap, A771, is associated with photon 
exchange, while the other, Ar] 2 , is associated with W exchange. 

At first sight, the probability for soft rescattering in such a process is rather small and we 
would expect S 2 ~ 1. The reason is that the transverse momentum, kt, distribution of the 
exchanged photon is given by the logarithmic integral 

for which the dominant contribution comes from the low k\ region. In other words, the main 
contribution comes from the region of large impact parameters, b, where the opacity of the 
proton is small. However, the minimum value of \t\, is not negligibly small. Note that the 
momentum fraction x p = 1 — £ associated with the upper proton can be measured with suffi- 
cient accuracy 18 , even without the tagging of the forward protons, by summing the momentum 



18 



The CDF collaboration [67] have demonstrated that this method provides an accurate determination of £. 
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(a) (b) 




Figure 13: Exclusive T production via (a) photon exchange, and (b) via odderon exchange, 
fractions 

6 = si m} + kl eW/vG (22) 

of the outgoing W and the hadrons observed in the calorimeters. As long as the gap Ai] 2 
is large, the dominant contribution to the sum £ = comes from the decay products of 
the W boson. For example, for rjw = 2.3(— 2.3), we expect an £ distribution centred about 
£ ~ 0.1(0.001). If we take Ar/ 2 > 3, then the cross section da(W + gaps)/<iln£ is 0.1(1) pb and 
~ 0.87(0.45) for £ ~ 0.001(0.1) and r] W = 2.3(-2.3), respectively. 

In the first LHC data runs it may be difficult to measure the absolute value of the cross 
section with sufficient accuracy. Most probably the ratio (VV+gaps/W inclusive) will be mea- 
sured first. In this case, the inclusive W production process (in the same kinematic region, 
Fig. 12(b)) plays the role of the luminosity monitor. Note that the cross section for inclusive 
W production is much larger than that with rapidity gaps. The reason is that an inclusive 
W is produced directly by qq fusion, which is prohibited for gap events since the colour flow 
produced by the t-channel quarks populates the A7712 rapidity gaps. 

Of course, the survival factor S 2 measured in VF+gaps events is quite different from that 
for exclusive Higgs production, which comes from smaller values of b. Nevertheless this mea- 
surement is a useful check of the model for soft rescattering. A good way to study the low 
impact parameter, b, region is to observe Z boson production via WW fusion, see Fig. 12(c). 
Here, both gaps originate from H^-exchange, and the corresponding b region is similar to that 
for exclusive Higgs production. The expected Z+gaps cross section is of the order of 0.2 pb, 
and 5' 2 =0.3 for Arji^ > 3 and for quark jets with E T > 50 GeV [68]. 

How can the LHC probe the unintegrated gluon f g distribution? Recall that the exclusive 
cross section, (12), depends on f g evaluated in terms of the integrated gluon density g(x,Q 2 ) 
for Q 2 t ~ 4 GeV 2 and x ~ M H /^. Now for Q\ = 4 GeV 2 and x = 10~ 2 the MSTW [69] 
and CTEQ [70] global parton analyses indicate that xg lies in the interval (3.2, 3.8). This is a 
big uncertainty bearing in mind that the exclusive cross section depends on (xg) 4 . To reduce 
the uncertainty associated with f g we can measure exclusive T production at the LHC. The 
process is shown in Fig. 13(a). The cross section for 7p — > Tp is given in terms of the same 
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Figure 14: The Rj distribution of exclusive two- and three-jet production at the LHC [73]. 
Without smearing, exclusive two-jet production would be just a 5-function at Rj = 1. The 
distribution for three-jet production is shown for two choices of the rapidity interval, 5rj, con- 
taining the jets; these distributions are shown with and without smearing. Here, we have taken 
the highest E T jet to have E T > 50 GeV. To indicate the effect of jet smearing, we have assumed 
a Gaussian distribution with a typical resolution a = 0.6/V-&T in GeV. 

unintegrated gluon distribution f g that occurs in Fig. 8. There may be competition between 
production via photon exchange 19 Fig. 13(a), and via odderon exchange, see Fig. 13(b). A 
lowest-order calculation (e.g. [71]) indicates that the odderon process (b) may be comparable 
to the photon-initiated process (a). If the upper proton is tagged, it will be straightforward to 
separate the two mechanisms; larger p t is expected for the odderon process of Fig. 13(b). 

A more detailed probe of f g may come from studies of central exclusive jet production 
at the LHC. We have already noted that the computation of f g in terms of the integrated 
distribution g involves the presence of a Sudakov-like form factor T, which ensures the absence 
of gluon emission in the evolution from Q t up to the hard scale ~ M^/2; and that T provides 
the infrared stability of the Qt integral in (12) for sufficiently large Ma- Also recall from 
Fig. 11 that a measure of exclusive dijet production at the Tevatron, pp — > p + jj + p, was 
obtained [60] by plotting the cross section in terms of the variable Rjj = Mjj/M^, where Ma 

19 The first experimental data on exclusive J/ip production obtained by the CDF collaboration [61] are con- 
sistent with the theoretical expectations for photon exchange [71], see also [72]. 
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is the mass of the whole central system. However, the Rjj distribution is smeared out by 
QCD radiation, hadronization, the jet algorithm and other experimental effects [60, 73]. To 
weaken the smearing it was proposed in Ref. [73] to study the dijets in terms of the variable 
Rj = 2Et (cosh f]*) /Ma , where only the transverse energy and the rapidity rj of the jet 
with the largest Et enter. Here 77* = i]j Ct — yA, where tja is the rapidity of the central system. 
Clearly, the largest Et jet is less affected by the smearing. 

In Fig. 14 we show the Rj distribution of both exclusive two- and three-jet production 
expected at the LHC For three-jet production we show predictions for two choices of the 
rapidity interval Srj within which all three jets must lie. If we take the largest Et jet to have 
E T > 50 GeV at the LHC, we see that the cross section for exclusive three-jet production reaches 
a value of the order of 100 pb. Of course, if we enlarge the rapidity interval 5r] where we allow 
emission of the third jet, then da/dRj will increase, see Fig. 14. Indeed, the measurement of 
the exclusive two- and three-jet cross sections as a function of E T of the highest jet allows a 
detailed check of f g and the Sudakov factor T; with much more information coming from the 
observation of the 5r] dependence of three-jet production. Note that the background from the 
inclusive interaction of two soft Pomerons (so-called double-Pomeron-e! xchange) should be 
small for Rj > 0.5, and can be removed entirely by imposing an E T cut on the third jet, say 
Et > 5 GeV. 

Finally, how can the LHC probe the gap survival probability, S^ nh , to enhanced rescattering, 
which violates soft-hard factorization? It appears at first sight, that we may study the role 
of enhanced absorption by observing the W+2 gaps process shown in Fig. 12(a). To do this 
one may vary the transverse momentum of the accompanying quark jet, q, and the size of 
the rapidity gap Ar^ 2 - For lower transverse momentum of the quark jet we expect a stronger 
absorptive effect, that should decrease with increasing A 772, since the number of partons in the 
rest of the rapidity interval increases. Unfortunately, the photon-initiated process of Fig. 12(a) 
occurs at large impact parameter b where the probability of rescattering is small. Moreover it 
will challenging to observe the quark jet at low Et in the relevant rapidity interval (—5 < 77 < 
—3). Therefore we will discuss other processes which depend on rescattering on intermediate 
pa! rtons. 

The observations we have in mind are the measurements of ratio R of diffractive (one-gap) 
events for W (or T or dijet) production as compared to the number of events for the inclusive 
process (shown in Fig. 12(b) for W production). These processes are shown schematically in 
Fig. 15. In other words, R is the ratio of the process in diagram (c) to that in diagram (a). 
That is 

no. of (A + gap) events a diff (xjp, g, /x 2 ) 2 2 

no. of (inclusive A) events a™ l (x = px F , /1 2 ) ^ cik ^ over fe ' ^ 

where a mcl and a dlff are the parton densities known from the global analyses of inclusive and 
diffractive deep inelastic scattering data, respectively. The heavy central system A is either W 
or a pair of high E T jets or T or a Drell-Yan pair. For W or ji + ji~ pair production the 
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Figure 15: Schematic diagrams for (a) the 
inclusive production of a system A, (b) and 
(c) for the diffractive production of A with- 
out and with 'enhanced' soft rescattering 
on intermediate partons. The system A is 
taken to be either a W boson or an T or a 
pair of high Et jets. 
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Figure 16: The predictions of the ratio R of (23) 
for the production of a pair of high E T jets. 



parton densities a are quark distributions, whereas for dijet or T production they are mainly 
gluon densities. Thus measurements of the ratio R will probe the gap survival factor averaged 
over the impact parameter h. 

To demonstrate the possible size of the effect, a simple model was used to estimate S^ nh 
[48]. The results for the dijet case are shown by the dashed curves in Fig. 16 as a function of 
the rapidity i/a of the dijet system. The enhanced rescattering reduce the ratios and lead to 
steeper i/a distributions, as illustrated by the continuous curves. Perhaps the most informative 
probe of S^ nh is to observe the ratio R for dijet production in the region Et ~ 15 — 30 GeV. For 
example, for Et ~ 15 GeV we expect S% nh ~ 0.25, 0.4 and 0.8 at i/a = —2, and 2 respectively. 



4 Conclusions 



We have described a model of soft high-energy pp (and pp) interactions with multi-components 
in both the s- and t-channels. The parameters of the model were adjusted to reproduce all the 
available data in the CERN-ISR to Tevatron energy range, and predictions are made for the 
LHC energy a/5=14 TeV. The inclusion of absorptive effects are vital. Low-mass diffractive 
dissociation is described in terms of a multichannel eikonal, whereas high-mass dissociation is 
included using a full set of multi-Pomeron-Pomeron diagrams. The full set is necessary as an 
analysis of triple-Regge data finds that the triple-Pomeron coupling, g 3P , is three times larger 
than that coming from the original analyses, which did not include absorption. A novel feature 
of the model is that the Pomeron is described in terms of 3 interconnected components with 
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different k t , which allows the soft to hard Pomeron transition to be studied. The soft (small k t ) 
component is heavily screened, so that dependent observables, such as the total cross section 
and the parton multiplicity at low k t , show relatively little growth with y/s; <r tot is prediced to 
be about 90 mb at the LHC energy ^/s=14 TeV. On the other hand, the large k t (or QCD) 
component suffers relatively little screening and the dependent observables have larger growth 
with y^i. 

It was emphasized that central exclusive production of a heavy system A, that is, the process 
pp — > p + A + p, offers an excellent opportunity to study the Higgs sector at the LHC via its bb 
decay channel, which, otherwise, is very challenging to observe. The crucial observation is that 
there exists a selection rule which greatly suppresses the pp — > p + bb + p QCD background. 
To identify the exclusive Higgs signal it is necessary to instal detectors for the outgoing very 
forward protons at a great distance (about 400 m) from the pp interaction point. These will 
enable the Higgs mass to be measured via the missing mass method to an accuracy 0(1) GeV. 
However, we must pay a price for the exclusivity of the Higgs signal. Soft rescattering can 
populate the rapidity gaps. We described how to calculate the survival probability of the 
gaps to both eikonal and enhanced rescattering. The latter, which involves rescattering on 
intermediate partons, breaks the soft-hard factorization of the amplitude for the process. Care 
must be taken to evaluate its effect. A vital step is the model of soft interactions with Pomeron 
exchange with different k t values, which enabled the b <-> k t correlation to be taken into account. 
It was found that the suppression caused by enhanced rescattering is numerically not large for 
the exclusive production of a heavy mass system (like a Higgs boson) at the LHC. 

We noted two results of NLO analyses of HERA data. First, 'flat' small-x gluons are 
obtained at the relevant scales in the global fits of deep inelastic and related data. Second, small 
values of the saturation scale Q s (b) are obtained for impact parameters larger than about 0.6 fm. 
These observations imply that the gluon density, which gives rise to the enhanced absorptive 
correction, is not large, and that the black disc regime is not reached at the LHC 20 . The cross 
section predicted for a 120 GeV SM Higgs boson, before acceptance and bb identification cuts, 
is 2 — 3 fb at the LHC energy ^/s = 14 TeV; with a signal-to-background ratio of 0(1). The 
cross section and signal-to-background ratio may be up to an order-of-magnitude higher in 
some SUSY Higgs scenarios. 

We detailed exclusive processes that have already been observed at the Tevatron with rates 
in agreement with predictions of the same model as that used for the Higgs estimates. This 
gives powerful support for the addition of forward proton taggers to enhance the discovery and 
physics potential of the ATLAS and CMS detectors at the LHC. Finally, we described a range 
of processes that may be observed in the early runs of the LHC which can provide valuable 
checks of all aspects of the theoretical formalism for exclusive processes. 

20 Note that the black disc regime of [66] was obtained using leading order (LO) gluons which grow steeply 
with 1/x. This growth is simply an artefact of the absence of the 1/z singularity in the LO quark splitting 
function. When NLO contributions are included the deep inelastic data are described by flat (or even decreasing 
with 1/x) gluons at the relevant low scales, see the discussion in [16]. 
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